Preprint typeset in JHEP style - HYPER VERSION 



TAUP 2929/11 



June 17, 2011 



Geometric scaling behavior of the scattering amplitude for DIS 



^ with nuclei 

G; 

3 ; 



b t 



Andrey Kormilitzin a * Eugene Levin" 6 t and Sebastian Tapia 

i— i ° Department of Particle Physics, School of Physics and Astronomy, Tel Aviv University, Tel Aviv, 

*5 69978, Israel 

Q-T h * 

Departamento de Fisica, Centro de Estudios Subatomicos, Universidad Tecnica Federico Santa Maria, 

^H, and Centro Cientfico-Tecnologico de Valparaiso, Casilla 110-V, Valparaiso, Chile 

jG : 

r ^ Abstract: The main question, that we answer in this paper, is whether the initial condition can influence 
on the geometric scaling behavior of the amplitude for DIS at high energy. We re-write the non-linear 
Balitsky-Kovchegov equation in the form which is useful for treating the interaction with nuclei. Using 

CO the simplified BFKL kernel, we find the analytical solution to this equation with the initial condition 

^sO given by the McLerran-Venugopalan formula. This solution does not show the geometric scaling behavior 
of the amplitude deeply in the saturation region. On the other hand, the BFKL Pomeron calculus with 

i—t the initial condition at xa = I/tuRa given by the solution to Balitsky-Kovchegov equation, leads to the 

* * i 

t> geometric scaling behavior. The McLerran - Venugopalan formula is the natural initial condition for the 
^ Color Glass Condensate (CGC) approach. Therefore, our result gives a possibility to check experimentally 
^ which approach: CGC or BFKL Pomeron calculus, is more adequate. 



Keywords: Color Glass Condensate, gluon saturation, BFKL Pomeron, calculus, non-linear evolution, 



geometric scaling behavior 



PACS: 12.38-t, 12.38. Cy,l 2.38.Lg, 13.60.Hd, 24.85. +p, 25.30.Hm 



* Email 
f Email 
'Email 



andreykl@post.tau.ac.il. 
leving@post.tau.ac.il, eugeny.levin@usm.cl. 
trockut@gmail.com 



Contents 



|]. Introduction. [I] 

General solution to the simplified non-linear equation. § 

pTT| Simplified equation for scattering with nuclei Q 

^ Simplified BFKL kernel | 

r < 1 | 



2.2.2 r > 1 



Traveling wave solution and the geometric scaling behavior of the scattering amplitude. |B| 

|3|. Dipole-nucleus amplitude: solution for one critical line and violation of the geometric 
scaling behavior. |^ 

Solution for two critical lines |ll] 

D] Equation for Y < Y A 



y solutions y 

|5|. Impact parameter dependence of the scattering amplitude 15 
ffl. Conclusions 16 



1. Introduction. 

Geometric scaling behavior of the scattering amplitude is one of the most qualitative and well founded 
features of the high energy scattering in the framework of the CGC approach. In the saturation domain 
it follows from the general structure of the non-linear equation (see Ref. [1]) that governs the scattering 
amplitude in high density QCD [2-7]. In the vicinity of the saturation scale the geometric scaling behavior 
was derived [8] from the linear evolution [9,14]. The geometric scaling behavior means that the amplitude 
turns out to be a function of one dimensionless variable: r = r 2 Q 2 s (x; b) instead of three variables: x - the 
fraction of the energy carried by the parton interacting with the virtual photon in DIS, r- the typical size 
of the interacting dipole (r oc 1/Q where Q is the photon virtuality), and b is the impact parameter of the 
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scattering process. Actually, the geometric scaling behavior reflects that the only dimensional parameter , 
that governs the scattering process in the saturation domain, is the saturation scale Q s and r is just the 
only dimensionless variable that we can construct. Implicitly, such a behavior assumes that the amplitude 
in the saturation region does not depend on the initial condition for the non-linear equation. As we have 
mentioned this idea looks consistent with the geometric scaling behavior of the scattering amplitude in the 
perturbative QCD region near to the saturation scale [8]. 

However, it has been noticed in Ref. [10] that the geometric scaling behavior of the scattering amplitude 
cannot be correct for the DIS with nuclei if the McLerran-Venugopalan formula [4] is used as the initial 
condition for this process. This initial condition itself shows the saturation behavior on the scale which is 
the saturation scale at low energy (Q s (x = x&b)). Generally , this behavior could be different from the 
geometric scaling one as it is shown in Ref. [10]. 

In this paper we revisit this problem and we try to answer the following questions: (i) could the initial 
condition affect the behavior of the scattering amplitude at r S> 1; (h) can we trust the McLerran- 
Venugopalan formula deeply in the saturation region; and (hi) what initial condition we need to use to 
reproduce the geometric scaling behavior. 

We believe that all these questions are needed to be answered not only due to pure theoretical inter- 
est but also because the geometric scaling behavior is seen experimentally both for hadron and nucleus 
scattering (see Refs. [11]). 



2. General solution to the simplified non-linear equation. 



The nonlinear Balitsky-Kovchegov equation for the scattering amplitude of the dipole with size r has the 
following form [6,7]: 



ON (r, Y; b) C F a s 



dY 2n 



2 



J d 2 r' K (r; r') x {2A f r ', Y;b — ^ (r - f')\ (2.1) 



N(r,Y;bj - N (r f , Y; b - ^ (r - r ') j N (r - r ', Y; b - 'J } 



where Y = ln(l/x) is the rapidity of the incoming dipole; is the imaginary part of the scattering 
amplitude and b is the impact parameter of this scattering process, as is the QCD coupling, below we will 
use the following notation as = asN c /ir. The BFKL kernel K (r; r') has the following form 

KM = (r ~_£)2r" (2 ' 2) 
2.1 Simplified equation for scattering with nuclei 

Considering the interaction of the dipole with a nucleus we can simplify this equation. Indeed, neglecting 
the non-linear term in Eq. ( |2.lD we have the linear BFKL equation and the amplitude for dipole-nucleus 
scattering can be written as 
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K FKL = J d 2 b'N* FKL (r, Y-b-b') T A (b') = T A (b) J d 2 W N N (r,Y; 6") = N N (r, Y; t = 0) T A (b) 

where b — b' = 6" and T A (b) is the number of the nucleons at given impact parameter b in the nucleus. It 
can be calculated as the following integral 



T A (b) = / dzp(b,z) (2.4) 

J — oo 

where p is the density of the nucleons in a nucleus and z is the longitudinal coordinate of the nucleon. This 
number depends on b but it is of the order of A 1 / 3 . 

In Eq. ( |2.3j ) we assumed that \b — b'\ < Rn *C 6' ~ i?^. One can recognize that this is typical 
Glauber-type assumption which holds for interaction with nuclei if the energy is not so high that the radius 
dipole-nucleon interaction can be considered as much smaller than the nucleus radius. Nn (r, Y;t = 0) is 
the imaginary part of dipole-nucleon amplitude at momentum transfer t = 0, which proportional to dipole- 
nucleon total cross section. 



Let us solve Eq. (|2.l|) using Eq. ( |2.3|) as the first iteration. One can see that the second iteration leads 



to the following expression 

dN { ? (r, F; t = 0) a s f d 2 r' r 2 t fS RFKT . . xN 

QY " = zf / (f - f/)2 r n >< I 2 ^(b) N* FKL (r', Y; t = 0)) (2. 

- T4(&)iV™(r,Y;* = 0) - N^ FKL (r', Y;t = 0)) N§ FKL (r — r F; t = 0) 

Introducing iV^ = N A /T(b) one can see that Eq. ( J2.5[ ) can be viewed as the first iteration of the 
equation 

a ^'Y Y ^ a) - If / =< W (r' y, * - «)) M 

- iV^ (r, F; t = 0) - T A (6) iV^ (r', Y; t = 0)) iV^ (r - f, F; t = 0) } 



The non-linear equation in the form of Eq. (|2.6| ) was firstly proposed in Refs. [2, 3] and it is valid 
for energies at which the radius of dipole - nucleon interaction is much smaller than the nucleus radius. 
The energy dependence of radius of dipole - nucleon interaction unfortunately cannot be solved in the 



framework of the Balitsky-Kovchegov non-linear equation of Eq. (2.1) (see Ref. [12]) since it demands 
some additional input from the non-perturbative QCD. The high energy phenomenology based on the soft 
Pomeron approach leads to 

R 2 (dipole - nucleon) = R 2 + a' P (0) Y <C R 2 A (2.7) 
in wide range of rapidity F including the LHC energies for Rq = WGeV~ 2 and Op(0) < 0.2GeV~ 2 . 



2.2 Simplified BFKL kernel 



The BFKL kernel of Eq. (2.2) is rather complicated and the analytical solution with this kernel has not been 



found. In Ref. [13] it was suggested to simplify the kernel by taking into account only log contributions. 
From formal point of view this simplification means that we consider only leading twist contribution to the 
BFKL kernel. Notice that the kernel of Eq. Q2.2j ) includes all twists contributions. We are dealing with 
two kinds of logs in Eq. (|2.l| ), which corresponds two different kinematic regions: r < 1 and r > 1. 

2.2.1 r < 1 

In this kinematic region we can simplify K (r;r r ) in Eq. (2.1) in the following way [13], since r' 3> r and 

J d 2 r' K (r, r') -> vr r 2 ^ (2.8) 

Introducing (r, Y; t = 0) = N^f (r, Y; t = 0) /r 2 we obtain 

92n \ (r,y; 'r°\ = f(^ f (r,Y-t = 0) - r 2 K 2 QCD T A (b) (nf)\r,Y;t = 0)) (2.9) 



One can see that the simplified kernel of Eq. fl2.8|) sums In ( r A.q CD 
2.2.2 r > 1 

The main contribution in this kinematic region originates from the decay of the large size dipole into one 
small size dipole and one large size dipole. However, the size of the small dipole is still larger than 1/Q S . 
This observation can be translated in the following form of the kernel 

d 2 r'K(r,r>) -+ vr f |^{f (2 .10) 

Jl/Qi{Y,b) r Ji/Ql(Y,b) \r-r\ 

One can see that this kernel leads to the In (r 2 Q 2 )-contribution. Introducing a new function N e J f (r,Y;t = 0) 
f dr 2 N e J f (r, Y;t = 0) /r 2 one obtain the following equation 



d 2 N e J f (r, Y; t = 0) \( „. dN e J f (r,Y;t = 0) \ ~ eff . , 



The Mellin transform of the full BFKL kernel of Eq. fl2.2|) has the form 



X(7) = / ^e-^if (r;r') = 2^(1) - tf>( 7 ) - </>(!- 7) (2-12) 
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where £ = ln(r jr ) and V'C- 2 ') = dlnr(2:)/cte with r(z) equal to Euler gamma function. The simplified 



kernel replaces Eq. (2.12) by the following expression 

for r > 1 



x(t) 



(2.13) 

for t < 1 ; 



1-7 

One can see that the advantage of the simplified kernel of Eq. ( |2.13[) is that it provides a matching with the 
DGLAP evolution equation [14] in Double Log Approximation (DLA) for r < 1. We will show below that 
this kernel leads to the geometric scaling behavior of the scattering amplitude. The other attempt [15] to 
use a simplified kernel is related to the BFKL kernel in the diffusion approximation, namely, 

1\ 2 ,Jf l x3N 



X ( 7 ) = u +D (j--J + 0M 7 --j I (2.14) 
with 

w = 4 a s In 2; D = 14a s C(3) (2.15) 

In this approach we loose any matching with the GLAP evolution. Both simplified kernels reproduce the 
geometric scaling behavior giving the illustrations to the general conclusions of Ref. [1]. 

2.3 Traveling wave solution and the geometric scaling behavior of the scattering amplitude. 

It is well known (see Refs. [1,2,15,16]) that the equation for the saturation scale does not depend on the 
particular form of the non-linear term in Eq. ( |2.1j ) and it has the form 

\n{Ql{Y)/Q\Y = Y )) = ^(A ( Y - y ) (2.16) 

with the critical anomalous dimension 7 cr given by 

O^fcr 1 - Tcr 

Inserting Eq. ( |2.12| ) in Eq. ( p.!7| ) one obtains 

lev = \ and ln(Ql(Y)/Q 2 (Y = Y )) = Aa s (Y-Y ) (2.18) 

In the vicinity of the saturation scale r — > 1 the behavior of the dipole amplitude has the form [8, 16] 

N(Y;r) oc (r^) 1 " 7 " (2.19) 

We illustrate this behavior approaching to the saturation scale from the perturbative QCD region (t < 1). 
In this region we can use Eq. (|2.9|) neglecting the non-linear term. This equation has a simple DLA solution 



eff 



exp (Y - Y ) In (r 2 Q* (Y = Y )) (2.20) 
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which leads to the following expression for dipole -nucleus amplitude (see Eq. fl2~3[ 



N A = T A (b) exp 2y/a s (Y - Y ) In (1/ p Q\ (Y = Y ))) + In (r 2 Q 2 (y = Y )) 



Therefore, in vicinity of the saturation scale if In r <C \/81n (Q 2 (y) /Q 2 (y = Yq)) the dipole-nucleus 
amplitude can be written as 

N A = T A {b)ct> (r 2 Q 2 ) 1/2 = foe** (2.22) 

where 

z = Aa s (Y-Y ) + In (r 2 Q 2 s {A; Y = Y )) = £ s + £ (2.23) 

where 

& = In (Q 2 (A; Y) /Q 2 S (A; Y = y ) ) = 4a 5 (y - y ) , and £ = In (r 2 Q 2 (A; y = y ) ) (2.24) 

The saturation scale for nucleus we defined as Q 2 S (A, Y = Yq) = T\ (6) Q s (N; Y = Yq) where Q s (N; Y = Yq 
is the saturation scale for the nucleon at the initial energy. <pQ is a constant that absorbers all pre-exponential 
factors in the DLA solution. It is instructive to notice that 4>q oc as [2, 13]. 

Inside the saturation region we are looking for the solution of Eq. Q2.11| ) in the form [10, 13]: 

N e J f = Tz\b)Jdg(l - e-^'< y )) (2.25) 

From Eq. Q2.25Q one can see that we can easily to calculate the dipole-nucleus amplitude 

N A = T A (b)N e N ff = 1 - e"^' y ) (2.26) 
Substituting Eq. ( |2.25| ) into Eq. ( |2.11| ) we obtain 

foe'* = a s N e N ff e^ (2.27) 
Canceling and differentiating with respect to £ we obtain the equation in the form: 

d 2 6 



dYd£ 

Using variable £ s and £ we can rewrite Eq. ( p. 27 ) in the form 

d 2 ct) 1 
dtsdt ~ 4 

or in the form of 



a s ( 1 - e"^) ) (2.28) 



1 - e~^ Y ^ ) (2.29) 
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for z defined in Eq. (|2.23| ) and x = £ s — £. 

Eq. ( 2.3C ) has general traveling wave solution (see Ref. [17] formula 3.4.1) 



J (j>o 



dtf/ 



c + 



KX + A z 



1 + e- 



(2.31) 



where c, (fro, A and n are arbitrary constants that should be found from the initial and boundary conditions. 



From the matching with the perturbative QCD region (see Eq. (2.22)) we have the following initial 
conditions: 

cj)(t = z = 0,x) = 0o ; 4>' z (t = z = 0,x) = \<j)Q (2.32) 



These conditions allow us to find that k = and c = for (fro <C 1. Therefore, solution of Eq. ( p. 31 ) 
leads to the geometric scaling since it depends only on one variable: z. It has the form [13] 



V2 



For arbitrary (fro the solution has the form 



(2.33) 



(2.34) 



saturation region 



perturbative QCD 



Solution of Eq. (2.33) we can use in the case of interac- 
tion of the dipole with rather dilute target. In this case the 
initial conditions are determined by perturbative QCD and 



have the form of Eq. (|2.32[) . 

One can see that if (fro oc a$ <C 1 the solution of 
Eq. ( P03p gives 

(fr = (fr e^ z (2.35) 



while at z » 1 Eq. fl2.33|) leads to 



I 

Figure 1: Saturation region: z = £ s +£ and x = 
£ s — £ for interaction with dilute target (proton) 



(2.36) 
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3. Dipole-nucleus amplitude: solution for one critical line and violation of the geometric 
scaling behavior. 



The main ingredient of Color Glass Condensate (CGC) approach is the assumption that there exists such 
value of energy that we can describe dipole-nucleus amplitude using the McLerran-Venugopalan formula [4] : 



N A (r 2 ;Y;b) = 1 - exp (- a 2 s Const r 2 ln (r 2 Q 2 s (Y = Y )) ) = 1 



exp 



r z Qi(A;Y = Y )) (3.1 



The last equation is a simplification of the original formula but it reflects the main physics of saturation 
and considerable simplify calculations. 

Eq. Q3.1| ) can be translated into the boundary conditions for <f) on the line Y = Yq (£ s = 0, see Fig. [|) 
that has the following form: 

<M£ S = 0;£) = ^ (3.2) 

while solution of Eq. ( 2.33 ) gives quite a different function at £ s = (see Fig. ^). Therefore we need to 
find a more general solution than it is given by Eq. ( |2.31| ). 

One of the general features of solution of Eq. ( 2.33] ) is 
the increase of <f> in the saturation region, it means that 
only in the vicinity of the critical line we have to keep term 
exp (—</>). Inside of the saturation region we can neglect this 
term reducing the equation to the simple one, namely, 



m 

5 - 




1 

4' ° r 



1 

4 



(3.3) 



Figure 2: Boundary conditions at Y=Yq. 
The red(upper) line is <f>(£) from McLerran- 
Venugopalan formula while the blue (lower) one 
is the solution of Eq. ( 2.33 ). 

boundary conditions and it has the form 

01 (z) 



d 2 (f> d 2 4> 
dt 2 dx 2 

with the initial and boundary conditions of Eq. ( |2.32| ) and 
Eq. (|3.2| ), respectively. 

It is well known that the solution of this equation is 
different for t = z < x (£ < 0) and t = z > x (£ > 0) [17]. 
For t = z < x (£ < 0) the solution is not affected by the 



-X + 



One can see that the general solution to Eq. (p.3|) has the form: 



(3.4) 



0(6,6 = + Fi&) + F 2 (0 (3.5) 

and the solution of Eq. ( p.3p can be obtained from Eq. ( |3.5| ) using the restriction from Eq. (2.32). For 
t = z > x (£ > 0) we need to take into account the boundary condition of Eq. (3.2). Using the general 
solution in the form of Eq. ( |3.5D and the matching condition on the line £ = 
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saturation region 



no GS / 



V 




Figure 3: Saturation region: z = £ s + £ and a; = £ s — £. 



0i (e = o) = 2 (c = o) 

simultaneously with the boundary conditions that has the form 



h (6 = o) 



we obtain the following solution for £ > 



62(z,0 = ^ 2 /8 - £ 2 /8 + o e € + 6 



00 



Therefore, the solution to the simplified Eq. ( [3. 3D has the following form 

01 (^) for £< 0; 

02 (2,0 for £ > 0; 
For the solution of the general Eq. (2.29) we have 

0(2; Eq. ( HI )) for £ < 0; 

0(6, £;) = < 

^ z 2 /8 - £ 2 /8 + 0(6; Eg- (UD) - 00 + 00 e« + i0o6 for e > 0: 
In Eq. fl3.9|) we assumed that for £ > we are approaching the solution of Eq. (|3. 3[) . 



(3.6) 



(3.7) 



(3. 



(3.9) 



(3.10) 
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One can see that solution of Eq. ( 0.10 ) does not 



show the geometric scaling behavior and solution of Eq. (2.29) 
depends both on z and £. It happens so due to the in- 
fluence of the boundary conditions. 

In Fig. ||]we plot the numerical solution* of Eq. ( 2.29j ) 
in the region £ > (see Fig. with the following 
boundary conditions: 

<f>(v = t„t = 0) = 0(6); 

0(y = & = 0,£) = 0o e*; (3.11) 

One can see that this solution does not show the geo- 
metric scaling behavior inside the saturation domain. 

To preserve the geometric scaling behavior we need 
to assume that for £ > at Y = Yq the behavior of the 
scattering amplitude is not given by the Glauber 
(McLerran-Venugopalan) formula but rather is deter- 
mined by the solution of Eq. ( |2.33D . 

The initial condition based on McLerran-Venugopalan 
formula stems from the main assumption that there ex- 
ists the rather low energy at which dipole rescatters 
in the nucleus but the emission of gluons can be ne- 
glected. At first sight, we do have arguments why the 
emission is small. Indeed, if a^A 1 ^ ~ 1 for a$Y < 1 (or 
Y <C 1/as ~ 1/A 1 ^) the emission of the gluon will be 
suppressed while the Glauber-type rescattering will be 
essential since the interaction with the nucleus will be 
proportional to as A 1 / 3 ~ 1. Therefore, we can choose 
the energy (x) which is large enough to use only the 
exchange of gluon for the dipole amplitude while the 
emission of gluon will be still suppressed. For this kinematic region we showed that the geometric scaling 
behavior of the amplitude is not valid. Since the McLerran-Venugopalan formula follows from the CGC 
approach and represents its key feature, we can claim that the CGC leads to the violation of the geometric 
scaling behavior in the kinematic region where r 2 Q s (Y = Yq) > 1. 

However high density QCD has two facets at the moment: the color glass condensate CGC approach 
[4-6,18] and the BFKL Pomeron calculus [2,3,9,19,20]. Both these approaches lead to the same non-linear 
Balitsky-Kovchegov equation [6, 7] for the dilute-dense system scattering in the large N c approximation 
which we consider here. However, in the BFKL Pomeron calculus the emission of gluons is taken into 
account even at small values of energy. In this approach the natural initial condition is Ntf (r, Y = Yq) = 



1.6 



1.4 



1.2 



0.8 



0.6 



0.4 



0.2 



Fig. |a 


' 4>(z£) 


S = 2 








i;<0 


) 0.5 1 1.5 2 


2.5 3 


z 




Fig. §b 





Figure 4: The exact solution of Eq. (2.2S) for func- 
tion 4> (z,y)(Fig. ||-a) and the dependence on z at 
different values of £ (Fig. [|-b). y — £ s , £ s = z — £ and 
0o is taken to be equal to 0.1. 



It is worth to mention that Eq. (2.29) has the form which does not depend on extra parameters and, using the numerical 



solution, we do not loose the generality of our approach. 
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r 2 Q S (N; Yq) (compare with Eq. ([3.1] )) and the value of Yq is much smaller that Ya = (1/3) In A. This case 
we consider in the next section. 

Before doing this we would like to draw your attention to the fact that the condition a^A 1 ^ 3 ~ 1 can 
be reached in QCD only for scattering of states with typical extremely short distances (say, onium which 
made of two very heavy quarks). Running QCD coupling for such states could be small of the order of 
a s r^> \/A l l®. For nuclei the typical as is determined by the size of nucleons and could be as small as 
as = 0.2-^0.3 but not smaller. In this case the situation changes crucially: summing all powers of asY will 
lead us to the BFKL contribution namely a s e w ° . This contribution can be larger than the re-scattering in 
the classical gluon fields. It happens so at large Ya = In A 1 / 3 since ujq is larger that 1/3 for as = 0.2 0.3 
(see Eq. (pT5|)). 



4. Solution for two critical lines 



4.1 Equation for Y < Ya 



this observation into Eq. ( |2.6| ) by changing the definition of Ta (b) in Eq. (2.4), namely, 

r+l/mx 



/ + L/171X 
dzp(b,z) 
- 1 /mi 



In the framework of the BFKL Pomeron calculus the rescattering with large rapidities but smaller than 
Ya = In A 1 / 3 should be treated using the non linear equation. In this kinematic region each dipole interacts 
with the number of nucleons that are smaller than pRa ~ A 1 / 3 and which actually is equal to e Y p/m 
where p is the density of nucleons in the nucleus and m is the proton mass [22,23]. We can incorporate 



(4.1) 



For i C xa = e Ya Eq. ( fP| ) reduces to Eq. ( p.4] ) while for 1 3> x S> xa Eq. (4.1) leads to e Y p/m for 
cylindrical nuclei. 

Introducing 

N A (Y,0 = T A (b;Y) N e J f (Y,0 
in stead of Eq. ( |2.26 ) and using Ta (b; Y) in the form: 

p 2R A oc A 1 / 3 for Y < Y A ; 



(4.2) 



T A (b;Y) 



(4.3) 



e Y p/m for Y < Ya; 



we can re-write Eq. (^^) in the following form: 



x < x A {Y > Y A ) : 



dN A (Y;i 



dY 



x > x A (Y < Y A ) 



dN A (Y;C) 



dY 



a S { / deN A (Y;^ - N 2 a (Y;^}; (4.4) 
N A (Y;i) = a s { / d?N A (Y;?) - N 2 a (y;^}; (4.5) 
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We will solve these two equations and these solutions should be matched on the line Y = Ya (see 
Fig. ||). These two equations have different critical lines. The critical line of the first one (see Eq. Q4.4D ) 
has been discussed in Eq. (2.16) and Eq. ( |2.17| ). It is shown as line 2 in Fig. [I] and has the form 



£ = In (r 2 Q 2 s (A-Y A ) 



£,2s 



4a s (Y - Y A ) 



(4.6) 



The easiest way to find the critical line for Eq. ( |4.5| ) is to search the solution to the general Eq. (|2.6| ) with 
T A (b) replaced by T A (b; Y) in the semi-classical form 



n a (y,0 



o cj(Y,OY-(l- 7 (Y;())t + S 



(4.7) 



This solution has a form of wave-package and the critical line is the specific trajectory for this wave- 
package which coincides with the its front line. In other words, it is the trajectory on which the phase 
velocity (v p h) for the wave-package is the same as the group velocity ( v gr ). The equation v gr = v p h has 
the following form fort Eq. 



Vph 



a s 



xiric 
i 



+ &s- 



1 



-asx' (7c 



■gr 



Solution to Eq. (4J3) gives 7, 



7cr 7cr 

^/as + 0{as) and it leads to the equation (see Fig. [[]) 
£ = -in = (1 + 2^)(Y A -Y) 



(4.8) 



(4.9) 




Figure 5: Two saturation regions for N: region I for Y <Y A = (1/3) In A and £ > £% s = (1 + 2^/ois)(Ya — Y); and 
region II Y > Y A = (1/3), In A and £ > &s = ~^a s (Y - Y A ). Here £ = In (r 2 Q S (Y = Y ; A)) where r is the dipole 
size 
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4.2 Solutions 



In both regions we will look for solutions using 



< 1 



(4.10) 



In region II we obtain Eq. ( |2.28 ) which we need to solve with the same initial condition as in Eq. ( gjH ). 
As we have discussed Eq. fl2.33|) gives the solution of this problem. 

In the region I using Eq. ( EL1Q1) we obtain after differentiation over £ 



1 - e" 



Na,yH ~ N M = ^sN A (1 - N' M j or y e-+ - 
Differentiating Eq. ( [4.11 ) over ^ we get 

4>ys, ~ 4>y - 4>i = -a s N A ^ + a s ( 1 - e 



a s N A e 



(4.11) 



(4.12) 



Eq. ( 4.1 2| ) has a simple solution for large Y and £. Indeed, assuming that is large in this region , 
Eq. ( [4.12 ) can be re- written in the form 



<PYzx ~ <Pz lZ1 ~ <PYP Z1 ~ 

witht z x = e + = (i + 2v^) (yk - y) + e 



y 21 



-OLSZ\(j) zl + a s 



(4.13) 



One can see that the common solution of the two following equations 

1. <?V 5 = a s ; 2. 0y + l = -a s zi; (4.14) 

will be the solution of Eq. ( 4.13j ). It is easily seen that such a solution has the general form 

(f>Ri = ast {Y A -Y) + a s (Y A - Y) 2 /2 + (Y -Y A ) + F (0 (4.15) 

where F (^) is the arbitrary function. The initial condition for this equation follows from the solution of 
Eq. (|2.22|) where Ta {b) is replaced by Ta (b; Y). They have the form 



(4.16) 



The boundary condition stems from the solution in region II (see Eq. (2.33)) and has the form 

4>m{Y = Y A ;i) = <f>(0 (4.17) 



Choosing F (£) = <j>{^) we see that solution 4>Ri matchers the boundary condition of Eq. ( 4.17 ) but 
not the initial condition of Eq. ( [4.16| ). We need to solve Eq. ( [4.12 ) in the region of small z\ to satisfy this 



'For the sake of simplicity we consider 1\f5ts <C 1 in this expression. 
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condition, but in this region we cannot neglect the term exp (—(j>) in Eq. ( P3 ). We can approach this 

(4.18) 



region solving Eq. (4.11) which can be rewritten in the form: 

0y (Y, z x ) - e^ Y ^ + 1 = a s N A (Y; Zl ) 
with 



N A (Y; Z1 ) = ] dz' x {l - e~< Y ^)) 
After differentiation of Eq. (4.18) with respect to z\ one obtains 

<h, zi - <j> Zl (Y, Zl ) e^ Y ^ = a s (l- e~^ Y ^) 



The initial and boundary conditions for Eq. (4.20) looks as follows: 



initial conditions: 
boundary conditions: 



<f> (Y, z x = 0) = <£ ; 
4>(Y = o, Zl ) = <f>U) ■ 



This equation has been solved numerically and the solution for <f> and Na = 1 — exp 
in Fig. 0. 



My, z) 




Fig. |-£ 



(4.19) 



(4.20) 



(4.21) 



is shown 



Figure 6: The exact solution of Eq. ( 4.20D for function Na {y,z) 
Cis = z i — C an d 4>o is taken to be equal to 0.1 and as = 0.2. 



Fig. g-b 

N(y,z) (Fig. g-a) and (Fig. @-b). y = £u, 



Therefore, the full solution has the geometric scaling behavior in the region II but shows the violation 
of the scaling behavior in region I as it follows from Eq. ( [4.15 ) and Fig. || 

Comparing this result with the conclusions of the previous section we see that the BFKL Pomeron 
calculus predicts the geometric scaling behavior in the saturation region for r 2 Q 2 S (A; Yq) > 1 and Y > Ya- 



- 14 - 



5. Impact parameter dependence of the scattering amplitude 

In this section we complete the study of the impact parameter dependence of the scattering amplitude 
that has been started in Ref. [24]. In Ref. [24] we claim that in the framework of the BK equation with 
the simplified kernel the impact parameter dependence can be absorbed into redefinition of the saturation 
scale, namely, 

For proton r = r 2 Q 2 (P;x) - 

or 

z proton 



'Q s (A;x) with Q s (A;x) = T\(b) Q 2 s {P;x) = t a for nuclei 

(5.1) 

^nucleus = ^proton + 21n R( 6 )) ( 5 - 2 ) 



Nil, b) 





5 6 7 

b in fm 



Fig. 0-b 



Figure 7: The exact solution of Eq. ( [2.34 ) for function N (z,b) = 1 — exp (— 4> (z, b)) (Fig. 0-a) for the interaction 
of the dipole with the gold and its dependence on b at different values of z (Fig. ^-b). y — £ s , £ s = z — £ and <fto is 
taken to be equal to 0.1. 



However, this claim is based on the solution of Eq. ( 2.33 ) which assumed that (j>o <C 1. In the general 
solution of Eq. ( [134]) one can see that 6-dependence cannot be reduced to changes in the value of the 
saturation scale. In Fig. [?] we plot the dependence of the scattering amplitude on the impact parameter 
using the realistic Ta (b) for the gold [25] and the value of (po in Eq. ( p. 22 ) taken from the fit of the HERA 
data [26]. The first glance at Fig. ^ shows that the typical value of b increases with z. It has a natural 
explanation. Indeed the width b Q of b distribution we can define as N (z,b = bo) = e -1 . Since the amplitude 
N has the geometric scaling behavior the value of bo can be determined from the equation 



N A (r A ) = N A M(bo) Q 2 s {P;x 



(5.3) 
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where Q s (P; x) is the saturation scale for the proton target. In other words, the typical 6 can be determined 
from the following equation: 

Ti(V) Q 2 s {P;x) = t a = Const (5.4) 



Recalling that 

one can see that 

where we use that 

and h is about 0.5 fm. 
In terms of z Eq. (| 



Q 2 s (P;x) = Q 2 o(^) X with A = 4a 5 



b = R A + -h\\n{l/x) 



„ , b>R A _ b -KA 
Ta{oo) > e h 



(5.5) 
(5.6) 
(5.7) 



has even a simpler form: 



h 

Ra + 7T z - 



2 ^proton 



(5.8) 



It is wort mentioning that Eq. (5.< 
saturation momentum. 



does not depend on the specific form of energy dependence of the 



It should be stressed that we obtain a logarithmic increase of the radius of interaction with z. It stems 
from Eq. ( |2.6| ) where we have integrated over the impact parameter of the nucleon. As we have discuss we 
can trust this 6-dependence in the kinematic region of Eq. fl2.7p . Certainly, for the dipole-gold scattering 
for z < 7 we can use this Glauber-type approximation. All problems with 6-dependence are originated from 
the large b dependence of the dipole-proton amplitude which falls down as 1/6 in perturbative QCD [12]. 
Implicitly we assumed that the non-perturbative corrections to the B-K equation has been taken into 
account for dipole- nucleon scattering in the transition from Eq. ( p.l[ ) to Eq. d2.6|) . In Fig. [8| we plotted the 
dependence on z the total cross section of dipole-nucleus intertaction for the case of gold 



Otot 



d 2 bN(z,b) 



(5.9) 



One can see a significant difference between realistic Wood-Saxon Ta (6) and the simplified one Ta (6) = 
p O (R A — 6) where p is the density of the nucleons in the nucleus. This picture demonstrates the significance 
of correct 6-dependence for calculation of the physical observables. 



6. Conclusions 

We hope that we answered three questions that we have discussed in the introduction. The first one: could 
the initial conditions affect the behavior of the scattering amplitude at r 3> 1. The answer is yes and 
we gave the explicit solution of Balitsky-Kovchegov equation which shows the violation of the geometric 
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Figure 8: Total cross section for dipole-gold 
amplitude a tot — 2 f d 2 b N (z,b) in the satu- 
ration region: the solution of Eq. ( 2.34 ) (curve 
1) and the solution to Eq. ( 2.29| ) with simpli- 
fied T A (b) oc 6 (R A - b) (curve 2). 



scaling behavior of the scattering amplitude if you use the McLerran - Venugopalan formula as the initial 
condition . 

The second question: can we trust the McLerran- Venugopalan formula deeply in the saturation region, 
has a kind of negative answer. In the sense that the McLerran- Venugopalan formula cannot be considered 
as the unique initial condition. We demonstrated that in the BFKL Pomeron calculus this formula should 
be replaced by the solution to the non linear equation in the region of Y < Y A = HA 1 / 3 ) for very heavy 
nuclei. This statement gives the answer to the third question: what initial condition we need to use to 
reproduce the geometric scaling behavior. 

It is well known that we have two approaches to high density QCD; the BFKL Pomeron calculus 
[2,3,19,20] and Color Glass Condensate [4-7]. Both lead to the same Balitsky-Kovchegov equation for 
DIS. The difference between them lays in the initial conditions. For the CGC the initial condition is the 
McLerran- Venugopalan formula which is valid in the classical gluon field approximation. On the other hand, 
for the BFKL Pomeron calculus a natural initial condition stems from the solution of B-K equation for 
Y < Ya- Therefore, we can formulate the main result of the paper in the following way. The CGC approach 
leads to the violation of the geometric scaling behavior for DIS with heavy nuclei for r 2 Q s (A; Y = Yq) > 1 
while the BFKL Pomeron calculus leads to the geometrical scaling behavior of the amplitude for Y > Y4. 
This result gives a possibility to check experimentally which of these two approaches is more adequate. 
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